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1  Institute  for  Nonlinear  Science 


I.  Introduction 


The  study  of  stratified  fluid  dynamics  has  primarily  concentrated  on  the  linear  or  nonlin¬ 
ear  interaction  of  familiar  modes  of  density  oscillations  in  the  internal  wave  frequency  range 
between  the  local  inertial  frequency,  /,  and  the  local  buoyancy  frequency,  N.  It  has  been  ob¬ 
served  by  Muller  and  his  coworkers,  Muller  (1988)  and  Muller  et  al.  (1986,1988),  that  in  the 
linearization  of  the  Eulerian  equations  for  internal  wave  motion,  in  addition  to  the  familiar 
internal  wave  oscillations,  there  is  also  a  zero  frequency  mode  which  is  clearly  identifiable  as 
geostrophic  motion  or  more  precisely  motion  which  possesses  potential  vorticity  and  whose 
space  and  time  scales  are  to  be  established.  In  these  papers  the  question  of  whether  there 
is  observational  or  computational  evidence,  Riley  et  al.  (1981),  for  these  potential  vorticity 
carrying  modes  has  been  addressed  at  some  length.  The  evidence  is  quite  strong  that  these 
modes  do  exist,  and  in  some  sense  of  equipartition  carry  as  much  energy  as  the  familiar,  zero 
potential  vorticity,  internal  wave  motions. 

The  purpose  of  this  paper  is  to  provide  the  consistent  formulation  of  the  dynamics  of 
these  modes  by  studying  their  interaction  in  some  detail  in  a  fully  Hamiltonian  context  using 
a  Lagrangian  specification  of  the  fluid  dynamics  of  a  stratified  fluid.  The  reasons  for  using 
the  Lagrangian  view  of  the  fluid  dynamics  are  twofold:  (1)  the  formulation  is  canonical  and 
all  the  familiar  tools  of  Hamiltonian  mechanics  can  be  employed  without  worry  about  the 
noncanonical  Poisson  brackets  which  arise  in  the  Eulerian  fluid  dynamics.  (2)  The  conserved 
quantities  which  appear  as  cyclic  coordinates  in  a  canonical  formulation  (in  particular  the 
potential  vorticity)  have  been  “reduced”  out  of  the  Eulerian  dynamics  by  construction,  and 
thus  they  appear  as  mysterious  or  at  least  obscure  consequences  of  the  Eulerian  equations 
of  motion.  As  cyclic  coordinates  in  the  Lagrangian  formulation,  one  may  treat  them  on  an 
equal  footing  with  the  interacting  internal  wave  modes.  Once  this  is  done,  the  equipartition 
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of  energy  among  these  modes  appears  as  a  familiar  consequence  of  similar  arguments  for 
equi-distribution  of  energy  among  modes  in  general  Hamiltonian  systems. 

A  word  about  “reduction”  is  in  order  here,  (Abarbanel  et  al.  (19S6)).  This  is  the  process 
whereby  one  goes  from  the  six  fields  in  a  Lagrangian  picture:  the  Lagrangian  positions 
and  their  conjugate  momenta,  to  the  four  fields  in  an  Eulerian  framework:  the  three  fluid 
velocities  and  the  fluid  density.  The  remaining  quantities  are  constants  of  the  motion,  and 
the  reduced  description  provided  by  the  Eulerian  formalism  restricts  evolution  of  the  fluid 
to  surfaces  where  these  quantities  are  constant.  Thus  the  Eulerian  formulation  is  unable 
to  shed  any  dynamical  light  on  the  evolution  of  these  degrees  of  freedom  and  leaves  us  in 
the  dark  about  how  the  degrees  of  freedom  in,  for  example,  potential  vorticity  interact  with 
the  internal  wave  modes.  The  Lagrangian  formulation  which  keeps  all  degrees  of  freedom 
overcomes  these  complications. 

Having  identified  the  modes  which  carry  potential  vorticity,  we  associate  them  with  both 
small  scale  oscillations  which  would  contribute  to  carrying  energy  in  the  internal  wave  fre¬ 
quency  range,  /  <  u  <  N  and  with  the  large  scale  mesoscale  motions  which  drive  the 
internal  waves  by  transfers  of  energy  to  them,  Muller  (1976, 19SS),  Muller  et  al.  (19S6,19SS) 
and  Brown  and  Evans  (1981).  [We  refer  to  these  modes  as  PV  modes  to  indicate  that  they 
carry  potential  vorticity  regardless  of  their  spatial  scales.]  This  formulation  then  allows 
us  to  investigate  the  interesting  question  of  the  transfers  of  energy  and  momentum  from 
mesoscale  geostrophic  modes  to  the  internal  wave  degrees  of  freedom  using  familiar  Hamilto¬ 
nian  mechanics.  The  problem  is  certainly  harder  to  formulate  in  Eulerian  terms.  If  the  initial 
conditions  of  the  problem  were  chosen  to  have  the  potential  vorticity  modes  concentrated  at 
small  scales  only,  then  the  interaction  among  the  modes  would  be  essentially  a  small  scale 
energy  transfer  is-ue,  though  eventually  the  energy  would  affect  the  larger  scale  motions  as 
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a  matter  of  principle.  The  flow  from  the  larger  scales  to  the  smaller  scale  internal  wave  and 
‘vortical’  modes  presumably  occurs  on  a  much  more  rapid  time  scale  than  energy  flow  in  the 
other  direction,  if  indeed,  the  latter  can  occur  at  all  before  being  dissipated  at  the  smallest 
scales  by  viscosity. 

We  also  formulate  the  problem  in  terms  of  a  compressible  fluid  and  would  recover  the 
incompressible  case  of  more  direct  physical  interest  by  taking  the  physical  limit  of  sound 
speed  going  to  infinity  at  the  end  of  any  calculations.  This  requires  us  to  carry  acoustic 
modes  of  the  fluid  motion  along  with  internal  waves  modes,  but  it  also  allows  us  to  ignore 
the  constraint  of  incompressibility  throughout  the  work.  Since  this  constraint  is  peculiar 
and  difficult  to  implement  in  the  Hamiltonian  framework,  it  is  convenient  to  have  it  out 
of  the  way.  Lighthill  (1978)  shows  in  a  linear  context  that  as  the  sound  speed  becomes 
large  compared  to  fluid  motion  speeds,  the  acoustic  and  other  modes  decouple,  as  one  would 
certainly  expect  on  physical  grounds. 

Our  work  proceeds  from  the  general  Lagrangian  formulation  of  the  fluid  dynamics  for 
three  dimensional  flows.  We  linearize  the  canonical  degrees  of  freedom  around  a  base  state 
of  no  flow  and  identify  the  quadratic  terms  in  the  Hamiltonian  for  the  problem.  After  a 
discussion  of  the  dispersion  relation  of  the  linearized  problem,  we  exhibit  a  canonical  trans¬ 
formation  from  the  original  canonical  coordinates  into  a  set  which  explicitly  possesses  the 
conserved  potential  vorticity  as  one  of  the  canonical  momenta  whose  canonical  coordinate 
is  absent  from  the  Hamiltonian.  That  is,  potential  vorticity  (PV)  appears  as  a  cyclic  coor¬ 
dinate  in  the  quadratic  Hamiltonian.  Next,  we  impose  periodic  boundary  conditions,  and 
identify  appropriately  scaled  Fourier  coefficients  as  canonical  variables.  It  is  then  possible 
to  explicitly  decouple  the  vortical,  internal  wave  and  acoustic  modes  by  a  further  canonical 
transformation.  In  other  words  the  quadratic  Hamiltonian  is  expressed  as  a  sum  of  inde- 
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pendent  harmonic  oscillator  Hamiltonians.  Finally  we  formulate  the  problem  of  interaction 
among  the  modes  identified  in  the  linear  problem  and  end  with  a  discussion  of  future  work 
including  both  numerical  directions  and  a  possible  way  to  utilize  the  full,  nonlinear  potential 
vorticity  as  a  cyclic  canonical  momentum  in  a  further  canonical  set  of  variables  for  the  prob¬ 
lem.  This  paper  treats  the  linear  modes  only.  Our  work  on  the  nonlinear  interaction  among 
the  modes  identified  here  will  appear  in  future  publications.  In  our  final  section  we  provide 
both  a  summary  of  the  work  in  the  paper  and  some  suggestions  on  the  work  to  be  done 
in  the  nonlinear  problem.  In  particular  we  note  that  the  potential  vorticity  in  Lagrangian 
formulation  is  quadratic  in  the  modes  we  identify  in  the  linear  problem.  This  is  in  contrast 
to  the  Eulerian  formulation  where  the  potential  vorticity  is  an  infinite  series  in  the  velocity 
and  density  perturbations  to  the  base  state.  This  suggests  that  working  to  all  orders  in 
the  nonlinear  problem  may  well  be  simpler  in  Lagrangian  picture  especially  if  one  wants  to 
emphasize  the  matter  of  potential  vorticity  conservation. 

II.  Action  Principle  in  Lagrangian  Formulation 

Our  starting  point  is  the  Lagrangian  for  the  fluid  dynamics  of  a  compressible  fluid  with 
an  energy  density  t{p)  and  canonical  coordinate  Y(r,  t),  which  is  the  Lagrangian  position 
vector  of  a  fluid  particle  with  label  r.  This  label  is  also  the  initial  position  of  the  fluid  particle 
for  we  require  Y (r,0)  =  r.  We  want  the  motion  to  take  place  in  a  rotating  frame,  so  we  add 
the  “rotational  potential-1  R(Y,t)  whose  curl  with  respect  to  Y  is  twice  the  local  rotation 
of  the  earth.  In  this  paper  we  will  work  entirely  in  an  f-plane  framework,  thus  we  can  take 

R(Y,!)  =  |(-xU  +  yy,)  (l) 

for  rotation  about  the  z-axis  at  // 2. 
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III.  Hamiltonian  for  the  Z  nearized  Motion 


We  want  to  isolate  the  linear  terms  in  the  equations  of  me* ion  which  arise  from  this 
Hamiltonian  and  the  canonical  Poisson  brackets 

{K(r,  i),  n,(r',  <)}  =  SijS3(r  —  r');  (8) 

i,j  =  1,2,3.  For  this  purpose  we  want  the  quadratic  terms  in  the  Hamiltonian  in  the 
expansion  of  the  dynamical  variables  around  the  base  state  Y(r,  t)  —  r  and  dtY(r ,t)  =  0 
using 

Y(r,  t)  =  r  +  X(r,  t),  (9) 

and 

n(r,t)  =  ~-i  x  r  +  7r(r,t).  (10) 

The  only  tricky  part  of  the  expansion  procedure  is  in  capturing  the  quadratic  terms  in  the 
energy  density  t{p)  —  e(pQd( r)/d(Y)).  For  this  we  need  the  result  that  when  y  =  r  +  X(r) 

d(Y)  „  .  d(X„Xt)  .  0(.V„A'3)  d(X,,X,)  d(X) 

aw  " 1  +  v'x  +  +  +  (u) 

This  leads  to  the  Hamiltonian  correct  to  second  order  in  the  coordinates  and  to  all  orders 


ffa(X,  Jr)  =  Jd3r  ^  -  X,x,]  +  X\  +  ■  X)2  -  pojA'sVj.  •  X±  , 


where  the  two  vector  Xj.  =  (AT,^),  and  the  speed  of  sound  is  defined  by  pp  =  c2  evaluated 
at  the  base  state  p  =  /)0(r). 

We  may  rewrite  this  form  for  the  quadratic  Hamiltonian  H?  in  the  following  suggestive 
fashion: 

tf2(X,  *)  =  j d3r  +  |(Xjn  -  X,x,)  +  of  XI  +  ^(Ap)2  +  \poN'Xl  ,  (13) 
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where  the  buoyancy  frequency  N  is 


and  the  acoustic  pressure  is 


9  dp0  _ 

pQ  dr3  c2  ’ 


A p  =  -p0c2  V  •  X  +  gpQX3. 


(14) 


(15) 


This  form  of  the  Hamiltonian  shows  clearly  the  separate  contributions  of  acoustic  energy 
and  internal  wave  energy,  see  Lighthill  (1978). 

The  equations  of  motion  following  from  this  Hamiltonian  and  the  Poisson  brackets 
{A',(r),7r,(r')}  =  <5u<53(r  -  r')  are 

dttX±  =  fdtXxx  z-VxC,  (16) 

and 

N2c2 

dttX3  = - —  (V-X)-a3C,  (17) 

where  (  =  Ap/ p0. 

It  is  immediate  from  the  taking  the  curl  of  the  Xx  equation  that  the  conservation  law 

dt  \d,x2  -  a2x,  +  /Vx  •  xx]  =  o  (is) 


holds.  This  is  the  linearized  form  of  the  conservation  law  found  in  Abarbanel  and  Holm 
(1987)  for  pressure  relationships  of  the  form  we  have  here;  namely,  pressure  a  function  of  p 
and,  through  po(r3),  also  a  function  of  r3.  That  conservation  law  states  that  the  following 
quantity  is  conserved: 


fl(r,. 


3  (dKJ_ 

"  h  \  dr2 


dY 


8K_a 

dr2  dr\ 


dt 


+  *n(Y,0 


(19) 


so 


50(r,  <) 
dt 


=  0. 


(20) 


This  is  the  potential  vorticity  as  shown  in  the  paper  of  Abarbanel  and  Holm  (1987). 

If  one  considers  the  buoyancy  frequency  constant,  then  these  linearized  equations  admit 
a  plane  wave  solution  of  the  form 


X(r,  t )  =  XP  exp|i(k  •  r  —  ut)], 


(21) 


with  X°  a  constant  vector.  The  frequency  is  related  to  the  wave  vector  by  the  dispersion 
relation 


jVV 


(a;2  -  /2)(^2  -  ik3(g  +  — ))  -  a,2jk|2c2  +  c2(*2/2  +  |*X|2.V2)  =  0. 
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This  has  the  approximate  solutions,  for  large  c,  of  an  acoustic  mode 


(22) 


and  of  internal  wave  modes 


uiw  — 


U'A  =  (kjc. 


{klf  +  \k±\2N2) 

|k|2 


(23) 


(24) 


There  is  also  a  zero  frequency  mode  associated  with  the  conservation  law,  and  in  this  mode 
we  have  the  condition  on  the  dynamical  variables 


JdtXj,  x  z  =  V±C,  (25) 

which  is  recognized  as  the  geostrophic  relation. 

One  of  our  goals  will  be  to  bring  out  the  conservation  law  as  resulting  from  a  cyclic  or 
ignorable  coordinate  in  the  Hamiltonian  by  making  an  appropriate  canonical  transformation. 
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The  canonical  transformation  which  isolates  the  conserved  quantity  as  a  canonical  mo¬ 


mentum  with  an  ignorable  coordinate  is  the  following 

$i(r)0  =  d\X2  —  d?Xu  (26) 

&(r,t)  =  Vl-X±,  (27) 

Q3(r,t)  =  A'3,  (28) 

Pi(r,t)  =  -  ftr,  +  ^[V±  •  Xa]),  (29) 

A(r,<)  =  ^4(Vi  ■  *j.  +  ^[a,x2  -  a,.v,]),  (30) 

P3(r,t)  =  tt3.  (31) 

Using  this  transformation  we  find  the  new  Hamiltonian 

Q,P)  =  [d3r  f^-{|VJ.P,|J+|VJ.P2|J  +  />!)-  /(?,/>,  (32) 

J  L^o 


+  +  ^~7f~{Q2  +  d3Q3)2  -  pogQiQs  ■ 

This  is  clearly  independent  of  Q j,  so  the  canonical  momentum  Px  is  conserved  <?(Fi(r,  t)  =  0. 
The  coupling  between  P\  and  Q 2  acts  as  a  constant  displacement  on  the  Q 2  motion.  We 
could  remove  it  by  a  noncanonical  transformation,  but  it  is  harmless  as  it  stands,  so  we  leave 
it  present.  Px  is  directly  related  to  the  conserved  quantity  one  deduces  from  the  equations 
of  motion  since 

Pi( r,  t)  =  M2  -  d2Xr  +  /V x  •  X±j  .  (33) 

Vi  1 

IV.  A  Modified  Dispersion  Relation 

We  now  exhibit  a  further  canonical  transformation  that  explicitly  decouples  the  internal 
wave,  acoustic  and  potential  vorticity  (PV)  motions  in  the  Hamiltonian  for  the  linearized 
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motion.  This  provides  a  physically  motivated  set  of  coordinates  in  which  to  study  the  effects 
of  the  nonlinear  interactions.  We  have  already  isolated  the  PV  degrees  of  freedom  from  the 
other  degrees  of  freedom  by  a  canonical  transformation  given  in  the  previous  section.  We 
shall  see  that  a  subtlety  is  involved  in  the  use  of  this  transformation.  This  will  be  discussed 
below  and  later  another  canonical  transformation  will  be  given  that  explicitly  decouples  the 
PV  and  internal  wave  modes  for  the  linearized  motion. 

In  order  to  accomplish  this  transformation  to  normal  mode  coordinates  it  is  necessary 
first  to  make  an  intermediate  change  of  variables.  This  is  related  to  a  peculiarity  of  the 
dispersion  relation  (22),  which  the  reader  may  have  noticed,  namely  the  occurrence  of  imag¬ 
inary  quantities.  This  might  suggest  an  instability  of  the  steady  state,  but  as  hinted  at  in 
Lighthill  (1978),  it  is  the  result,  of  not  looking  at  the  problem  in  quite  the  right  coordinates. 
A  transformation  to  a  set  of  proper  coordinates  will  give  a  well  behaved  dispersion  rela¬ 
tion,  with  four  real  frequencies  at  each  choice  of  wave  number.  These  correspond  to  pairs 
of  acoustic  and  internal  waves  travelling  in  opposite  directions.  The  form  of  the  Hamilto¬ 
nian  in  Eq.  (12)  itself  suggests  the  correct  transformation  by  the  following  reasoning.  One 
would  like  to  make  the  problem  of  finding  normal  mode  coordinates  algebraic  by  using  a 
Fourier  representation  of  the  degrees  of  freedom.  The  appearance  of  factors  of  Po(r)  prevent 
a  straightforward  implementation  of  this  plan.  These  factors  go  away  by  using  the  following 
canonical  transformation 

v(r,  0  =  /5<T3/2(r)7r(r,  <),  x(r,  t)  =  p01/2(r)X(r,  t).  (34) 

In  terms  of  these  variables  the  quadratic  Hamiltonian  is  expressed  as 

•^2(x,v)  =  J d3r  ^v2  +  ^(x2v1  -  XiV2)  +  (35) 

C?  a/*^ 

+  ^^)2  +  (~-9)^x-^  +  ^-x2  , 
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where 


_  N2  t  g  d\ogp0(r3) 

a= - h-r  = - j - •  (36) 

g  c2  dr3 

Now  the  quadratic  Hamiltonian  contains  no  explicit  po  factors  so  that  a  transformation  to 
Fourier  modes,  assuming  N  to  be  constant,  will  indeed  make  the  normal  mode  problem 
algebraic.  Before  we  go  on  to  this,  it  is  interesting  to  compute  the  dispersion  relation  again, 
using  Hamilton’s  equations  for  x  and  v.  We  easily  deduce 

Ot(? 

dtlxx  =  fdtxxxz  +  Vx  c2 V  •  x  +  (— - g)x3  ,  (37) 

dttxz  =  c2a3( V  •  x) - — z3  -(— -5)Vx-xx, 

=  -^fx3-(~-g)V-x  +  d3  c2V  •  x  +  (^  -  g)x3  .  (33) 

This  yields  the  dispersion  relation  for  plane  waves 

(w!  -  /2)(u,J  -  ^f))  -  u.3|kpc!  +  S(klp  +  Ifcj.pJV2)  =  0,  (30) 

where,  as  before,  there  is  also  a  zero  frequency  solution  corresponding  to  the  geostrophic 
motion.  Comparing  this  dispersion  relation  for  combined  acoustic  and  internal  wave  modes 
to  the  previous  dispersion  relation,  Eq.  (22),  we  see  that  the  following  modification  has  taken 
place 

ih  — ♦  j.  (40) 

We  shall  see  later  that  the  dispersion  relation  in  its  new  form  has  four  real  roots. 

V.  Normal  Mode  Coordinates 

We  now  take  our  fluid  to  be  confined  to  a  periodic  box  with  horizontal  sides  of  length  L 
and  height  D.  Our  coordinates  and  momenta  can  now  be  expressed  as  Fourier  series 

*.(M)  =  yj==  Hg.(n,f)exp(t(Anj.  •  r±  +  ftn3r3)),  (41) 
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Expressing  the  quadratic  Hamiltonian,  tf2(x,v)  Eq.  (35),  in  terms  of  these  Fourier  coef¬ 
ficients,  we  obtain 

Hi (q,p)  =  p(n)-p(-n)  +  /(g2(n)p1(n)-g,(n)p2(n))  (46) 

f 2 

+  +  c2(An-i-  •  9x(n)  +  «”393(n))(Anx  •  gx(-n)  +  Kn3q3(- n)) 

?  2  2  1 
QCZ  Cl  C 

+  *A(—  -5)(nx-gx(n)93(— n)-nx-gx(-n)93(n))  +  -^-g3(n)g3(-n)  . 

The  evolution  equations  for  the  Fourier  coefficients  are  Hamilton’s  equations 

dgi{ n)  dH2  dpi( n)  _  dH2  ,47* 

dt  dpi(  n)1  dt  8q,(n) 

It  would  now  seem  that  in  our  search  for  normal  mode  coordinates  we  should  first  im¬ 
plement  a  Fourier  version  of  the  canonical  transformation  in  Eq.’s  (26)-(31),  which  isolates 
the  PV  motion.  However  there  is  an  issue  here  which  we  mentioned  above.  An  examination 
of  the  canonical  transformation  reveals  that  it  is  not  well  defined  for  phase  space  variables 
independent  of  rx  because  of  the  occurrence  of  the  inverse  perpendicular  Laplacian  in  the 


transformation.  This  is  a  warning  that  in  the  Hamiltonian  (46)  we  should  treat  the  case  of 
mode  amplitudes  satisfying  ni  ^  0  separately  from  the  case  nx  =  0.  Indeed  the  isolation  of 
the  PV  motion  can  be  performed  by  a  Fourier  version  of  the  canonical  transformation  (26)- 
(31)  for  the  case  nj.  ^  0.  The  case  ni  =  0  must  be  treated  differently,  however  normal  mode 
coordinates  are  still  easy  to  find. 

In  accordance  with  these  observations,  we  write 

Hi{q,p)  =  tfi(g(n),p(n))  +  Hv(q{0,0,n3),p(0,Q,n3))  (48) 

where  H±  is  the  part  of  the  Hamiltonian  involving  a  sum  over  Fourier  coefficients  satisfying 
nx  ^  0,  and  Hv  (V  for  “vertical”)  is  the  remaining  part  of  the  Hamiltonian  with  nx  =  0. 

A.  The  Case  of  Vanishing  Horizontal  Wavenumber  (n_L  =  0) 

For  convenience  we  write 


g(0,0,n3) 

=  q{n3), 

(49) 

p(0,0,n3) 

=  P{n3). 

(50) 

Then  from  Eq.  (46),  we  may  write 


Hv{q,p)  =  111 

L  r>3 

P 


Px(n a)  •  Pi(“  ”3)  +  /(?2(n3)Pi(n3)  -  9i(n3)p2(n3)) 


+  -jq±M  ■  9i(~”3)  +  P3(n3)p3(-n3)  +  7Z3q3(n3)g3(~n3) 


(51) 


where 


7ns  =  «2n3C2  + 


2  2 
a  c 


(52) 


For  these  nx  =  0  modes  the  motion  in  q3,  p3  is  already  explicitly  decoupled  from  qx,  p±. 
It  is  clear  that  the  evolution  of  the  q3,  p3  represents  acoustic  waves  with  allowed  frequencies 
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u>nj  =  ±7„3.  A  canonical  transformation  that  diagonalizes  the  remainder  of  Hy  is 


apv(n3)  =  ^<7i(n3)  +  jp2{-n3) 

(53) 

bpv(n3)  =  -tq2(-n3)  +  pi(n3) 

(54) 

=  ^?2(~n3)  +  jPi{n3) 

(55) 

bi(n3)  =  -tqi{n3)  +  p2{-n3). 

(56) 

and  for  uniformity  of  notation  we  will  also  set 

“A(n  3)  =  <?3(n3) 

(57) 

^(n3)  =  p3(n3). 

(58) 

The  notation  PV,  I,  and  A  indicates  potential  vorticity,  internal  wave,  and  acoustic  modes 
respectively.  It  can  be  checked  easily  that  the  a’s  and  the  6’s  form  conjugate  pairs 

{a:(^3)j  5j'(m3)}  —  (59) 

with  the  remaining  brackets  vanishing.  The  indices  i  and  j  now  run  over  PV,  I,  and  A.  In 
terms  of  these  variables  we  have 

Hv{ a,b)  =  ^  ]QMn3)M-n3)  +  f2ai{n3)ai(-n3) 

+  ^(n3)^(-n3) +  7^a^(n3)a/1(-n3)].  (60) 

We  see  that  the  PV  degrees  of  freedom  are  completely  absent  from  Hv,  and  the  allowed 
internal  wave  frequencies  are  simply  ±/.  These  results  are  in  agreement  with  the  fact  that  / 
is  the  minimum  allowed  linear  internal  wave  frequency  and  is  reached  when  kx  =  Anx  =  0. 
PV  or  geostrophic  motions  vanish  when  nx  =  0.  We  now  proceed  to  the  more  complicated 
case  when  nx  ^  0. 
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B.  The  Case  of  Nonvanishing  Horizontal  Wavenumber  (nx  0) 


The  part  of  the  Hamiltonian  involving  a  sum  over  the  nx  ^  0  modes  is 


=  \Yl 

&  n 


f 2 

p(n)  •  p(— n)  +  f(q2(n)pl( n)  -  9l(n)p2(n))  +  ^-9x(n)  *  9x(~n)  (61) 


+  c2(An±  •  9x(n)  +  Kn3q3(n)){\n±  ■  q±{-n)  +  nn3q3(- n)) 

Cl(?  Cp  C2 

+  iX(—  ~  ^)(ni.  *  9x(n)?3(-n)  -  nx  •  <7x(— n)<?3(n))  +  -£-q3{n)q3(-n) 


The  prime  on  the  summation  symbol  indicates  that  the  terms  nx  =  0  are  omitted  from  the 
sum.  A  Fourier  version  of  the  canonical  transformation  in  Eq.’s  (26)-(31)  will  now  isolate 
the  PV  motion  from  the  internal  wave  and  acoustic  wave  motions.  This  Fourier  version  of 
the  transformation  may  be  written 


<2x(n) 

=  fA(n192(n)  -  n29i(n)), 

(62) 

<?2(n) 

=  iA^i^InJ  +  n^n)), 

(63) 

<?3(n) 

=  ?3(n), 

(64) 

^i(n) 

—i 

A  nx2 

«iP2(n)  ~n2Pi(n)  + 

^(«i9i(-n)  +  n292(-n)) 

(65) 

P2(n) 

—i 

A  nx2 

rciPi(n)  +  n2p2(n)  -f 

t(niq2(-n)  -  n2?1(-n)) 

(66) 

P3(n) 

=  Ps(n). 

(67) 

From  these  formulae  it  is  evident  that  the  nx  =  0  modes  did  require  separate  treatment.  It 
is  again  simple  to  verify  that  the  Q's  and  P's  form  canonically  conjugate  pairs 


{<2i(n),P;-(m)}  =  Mnm, 


(68) 


with  the  remaining  brackets  vanishing.  In  terms  of  these  coordinates  we  have 


+  Tm<Mn)CM-n)  +  i;,Q3(n)(?3(-n)  -  2/<J2(n)P,(n) 

A  ni* 

+  c(<Ti  +  i<T2(n3))Q3(n)Q2(-n)  +  c(ctx  -  t<T2(n3))Q2(n)Q3(-n)],  (69) 


where 


72 

/n3 


=  K7n\c '?  + 


a2c2 


x  =  /2  +  A2nx2 


qc  g  .  , 

o-i  =  — - ,  tf2(n3)  =  «n3c.  (<0) 


Once  our  intermediate  canonical  transformation  given  by  Eq.  (34)  is  accounted  for,  P) 

given  in  Eq.  (69)  is  the  Fourier  space  version  of  P2(Q,P)  given  in  Eq.  (32).  We  see  that 
Qi(n)  does  not  occur  in  Hx{ Q,P),  so  the  Pi(n,f)  are  conserved  for  each  n.  Therefore  Qx 
and  Px  describe  the  PV  degrees  of  freedom.  The  term  involving  Qj(n)Pi(n)  in  Eq.  (69) 
appears  to  couple  the  PV  motion  to  the  internal  wave  and  acoustic  modes.  The  constancy 
of  the  Pi(n)  means  that  this  term  represents  just  a  constant  displacement  of  either  the 
acoustic  or  internal  wave  modes.  These  displacements  can  be  removed  from  the  quadratic 
Hamiltonian  by  a  noncanonical  change  of  variables,  but  we  choose  to  carry  them  along,  see 
Arnol’d  (1978).  We  now  give  one  further  canonical  transformation  that  serves  to  explicitly 
decouple  the  internal  wave  and  acoustic  modes 


AA(n) 

Ai{n) 


1 

A ^ 


(<M»)  +  *("))  • 

+(n)  -  wl(n)  \  c{<rx  -  ia2{n3))  J 


72, -ul(n)  <4(n)  -  7n, 


1 

Anj.  \l  w2  (n)  —  w2  (n) 


Q2( n)  - 


c(ax  -  icr2(n3)) 


Qs( n)  , 


B  (n)  = 

A  \/(wl(n)  “  wl(n))(u.’2  (n)  -  723) 

x  ((«+(n)  ~  ln3)Pi{n)  +  c(<rx  -  icr2(n3))P3(n)), 
B,(n)  =  .  ,  Anj- 

X  ((7n,  -  wi(n))P2(n)  -  c(<TX  -  i'cr2(n3))P3(n)), 


(71) 


(72) 


(73) 


(74) 
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where 


a4(n)  =  pl  +  -ri,  ±  m± + ii?  -  4  + a-xl  -  jv2))]1'1 

=  PL  +  7nj  =  VWni  _  Ti)2  +  4A2C3Tlj.2|<7|2,  (75) 

and  where  a  =  ax  4-  i<r2.  Note  that  |<rj2  =  723  —  N2.  We  see  from  this  that  both  are  real 
and  positive.  Also  we  see  that  the  quantities  723  —  u2_  and  a;2  —  723  are  positive.  All  the 
square  roots  in  the  canonical  transformation  are  therefore  well  defined.  Again  for  uniformity 
of  notation  we  also  put 


A/>v(n)  =  Qi(n),  Bpy(n)  =  Px{n).  (76) 

This  transformation  may  appear  mysterious,  however,  it  can  be  motivated  by  a  standard 
normal  mode  analysis  of  the  linear  equations  of  motion  for  Q  and  P.  u>±  are  simply  the  solu¬ 
tions  of  the  quadratic  equation  for  u2  in  the  dispersion  relation  (39),  where  the  wave  numbers 
have  been  properly  “quantized”  through  the  imposition  of  periodic  boundary  conditions 


k  =  (An1?  An2,  nn3). 


The  A's  and  the  B's  form  canonically  conjugate  pairs 


{A,(n),  5j(m)}  —  <5,j<5nm, 


with  the  remaining  brackets  vanishing. 


Note  that  as  the  speed  of  sound  becomes  infinite, 


-  |k|V  =  (A2n2  +  K2n\)c\ 


and 

2  k\N2  +  k2f2 
k\  +  k2 


(77) 


(78) 


(79) 


(80) 
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Now  we  express  the  nj.  ^  0  part  of  the  Hamiltonian,  Hx  in  terms  of  the  normal  mode 


coordinates 


Hd  A,B) 


“ 


\2n±2Bpv(n)Bpv{—n)  +  n)  +  B/(n)£/(— n) 

+  (n)^(n)A/t(-n)  +  wl(n)/l;(n)yl/(-n) 

-  2f\n±Bpv(n) 


(81) 


^j(n)  -  72, 

^ul(n)-ui{n) 


^(n)  + 


72  —  u>2  (n) 


\  w2(n)  -u;l(n)‘ 


From  Hamilton’s  equations 

dAi(  n)  dHd  A,B)  d£,(  n)  dHx{  A,B) 

df  “  di 5i(n)  ’  dt  ~  8Ai(  n)  ’ 

we  can  now  deduce  the  normal  mode  equations  for  the  evolution  of  either  the  momenta 
B(n,t)  or  the  coordinates  A(n,f), 


ApV{n,t)  =  \ 2nx2BpV(n,t)  -  f\nx 


u>i(n)  —  7  2 

+  V“/  In  j a  f 

-  —AA{r\,t) 


V2(n)-w!(n)‘ 


(82) 


+ 


-"-(n) 


A/(n,t) 


\  w2(n)  -w2(n) 
Bpv{n,t)  =  0, 

AA(n,t)  =  B^(n,f), 

B\{  n,t)  =  -w2(n)^(n,<)  +  /Anx 

A;  =  £;(n,f), 

Bj(x\,t)  =  -w2  (n)/l/(n,f)  +  /Anx, 


o’ 2  (n)  —  72 

+■  -  -BPv(n,t), 


w2(n)  -u£(n)‘ 


7n2,~^(n) 

N  (n)  “  “»i(n) 


Bpv(n,t). 


So  we  do  indeed  have  normal  mode  coordinates. 
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VI.  Summary;  Future  Work 


The  work  in  this  paper  has  concentrated  solely  on  the  identification  of  the  linear  normal 
modes  for  the  evolution  of  an  inviscid,  stratified  fluid — all  done  in  Lagrangian  formulation. 
We  have  worked  with  a  compressible  fluid,  so  the  three  kinds  of  linear  modes  are  an  acoustic 
mode,  internal  wave  modes,  and  a  low  frequency  potential  vorticity  carrying  mode  which 
satisfies  the  geostrophic  balance.  The  Hamiltonian  formulation  of  the  problem  in  Lagrangian 
fluid  variables  provides  both  a  fully  diagonalized  Hamiltonian  for  the  linear  problem  and  a 
framework  for  the  interaction  of  these  linear  modes  in  a  fully  nonlinear,  fully  Hamiltonian 
fashion. 


To  provide  a  summary  in  one  place  of  the  developments  and  various  canonical  transforma¬ 


tions  developed  in  this  paper,  we  collect  the  results  here.  First  of  all  we  have  the  Lagrangian 
coordinates  Y(r,  t)  =  r  +  X(r,  t)  and  canonical  momenta  iT(r,  t)  =  zxr  f/2  +  7r(r,  t)  which 
satisfy  the  canonical  Poisson  bracket  relation  {X,(r),  7rJ(r')}  =  StJS3(r  —  r')  with  all  other 
brackets  vanishing.  The  Hamiltonian  in  these  canonical  coordinates  takes  the  form 


H(X,  tt)  =  J  d3r  ■ 


7T-p0R(X)|2 


+  pogX3  +  Po^(-~j)  , 


where  the  quantity  £  is 


,  d(xux,)  ,  d(XuX,)  .  d(X,,X3)  .  3(X) 

?  =  1  +  V'X  +  ’«(^  +  '3(^  +  ^(^  +  'W-  1‘ 

The  expansion  of  this  Hamiltonian  to  second  order  in  powers  of  X(r,  t)  was  found  to  be 


'<*•*>  =  /  *  [ 


+  {|AV,  -  JVjttj]  +  £2pxl  +  ■  X)J  -  poSX3Vj.  ■  Xj.  , 


using  dpo(r3)/dr3  —  —gpo^s)  and  po  =  pleP{po).  This  is  the  quadratic  Hamiltonian  we 
diagonalized  in  the  body  of  this  paper.  To  do  so  we  performed  various  canonical  transforma¬ 


tions  from  the  phase  space  coordinates  X  and  7 r  to  others  which  evolve  as  simple  harmonic 
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oscillators.  In  a  moment  we  will  make  the  complete  connection  between  these  normal  mode 
coordinates  and  the  X  and  7r  coordinates. 


We  may  use  the  normal  canonical  coordinates  in  the  full  Hamiltonian  H(X ,  7r)  written 
above  and,  with  Hamilton’s  equations,  have  a  set  of  fluid  evolution  equations  which  are  not 
amenable  to  analytic  solution  but  can  be  addressed  numerically  with  the  advantage  over 
the  formulation  in  X,  tt  form  that  we  know  exactly  the  linear  interpretation  of  each  of 
the  modes:  acoustic  (A),  internal  wave  (I),  and  potential  vorticity  (PV)  as  they  enter  the 
Hamiltonian.  The  PV  modes  enter  in  a  nice  way  in  the  quadratic  Hamiltonian:  only  the 
PV  momenta,  called  Bpv(r,t),  are  present.  Their  canonical  conjugates,  called  Apy(r,t),  are 
absent,  so  the  Bpv  are  constant  in  time.  The  Apv(r,  t)  are  present  in  the  potential  energy 
term  Pc9^3  +  This  means  that  the  Bpv(r,t)  will  not  be  conserved  in  the  nonlinear 

interactions.  It  is  interesting  to  note  that  the  conserved  momentum  B  >v(r,t)  is  found  in 
the  potential  energy,  so  the  dynamics  of  regions  of  the  fluid  with  different  potential  vorticity 
will  generally  be  different.  This  would  seem  consistent  with  some  of  the  observations  of 
Sommeria  et  al.  (19S9). 

One  may  regard  this  as  a  problem,  but  it  is  not  a  serious  one.  This  is  for  two  reasons. 
We  know  from  general  considerations,  Abarbanel  and  Holm  (1987),  the  exact  form  of  the 
conserved  potential  vorticity  fl(r ,f) 


Po(r3)n(r,<) 


k  \  Or,  dr2 


dYn  dUn  ) 
dr2  dri  J  ’ 


(86) 


and  this  is  proportional  to  Bpv(r,t)  in  linear  order.  The  full  f!(r,  t)  is  quadratic  in  the  coor¬ 
dinates  and  can  be  easily  tracked  throughout  any  investigation  of  the  nonlinear  interactions. 
In  other  words,  the  exact  potential  vorticity  can  be  identified  in  the  interaction  and  in  the 
solution  of  the  equations  of  motion. 
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There  is  another  attractive  possibility  which  we  will  investigate  in  our  subsequent  papers 
on  this  subject.  The  full  potential  vorticity  is  only  a  quadratic  polynomial  in  canonical 
coordinates  It  may  be  possible  to  make  a  canonical  transformation  in  which  fl(r)  becomes 
a  new  “momentum”  with  its  conjugate  “coordinate”  absent  from  the  Hamiltonian.  Then  the 
conservation  law  will  be  evident  but  the  presence  of  fl(r)  in  the  interaction  potential  energy 
will  still  provide  a  driving  of  internal  wave  and  acoustic  motions  by  PV  states. 

One  may  use  this  observation  to  set  up  an  interesting  problem  where  working  to  some 
order  in  the  nonlinearity  one  specifies  the  PV  modes  at  t  =  0  knowing  they  do  not  evolve 
(on  an  f-plane).  They  then  act  as  time  independent  driving  terms  for  the  other  modes.  If 
one  establishes  a  mesoscale  flow  at  an  initial  time,  its  energy  will  transfer  to  internal  wave 
motions  as  the  system  evolves.  The  rate  of  transfer  can  be  compared  to  the  rates  observed 
by  Brown  and  Evans  (1981). 

The  decoupling  of  the  acoustic  modes  when  c2  — ►  oo  is  certainly  expected  on  physical 
grounds,  but  we  have  not  exhibited  this  explicitly.  While  we  will  deal  with  this  in  our 
future  work,  we  have  two  observations  which  may  be  useful  in  this  regard.  First,  there 
actually  may  not  be  a  smooth  transition  in  a  formal  sense  to  infinite  sound  speed  since  the 
Hamiltonian  formulation  of  that  problem  ab  initio  requires  introducing  the  pressure  as  a 
Lagrange  multiplier  which  serves  to  enforce  the  constraint  that  the  Jacobian  d(Y(r,  t))/d(r) 
always  equal  unity.  This  departs  from  the  straight  line  Hamiltonian  path  we  have  developed 
here  and  requires  complications  to  introduce  this  constraint  in  Poisson  brackets  and  other 
necessary  formal  apparatus,  Abarbanel  et  al.  (1986).  If  one  has  to  avoid  the  limit  in  the 
Hamiltonian  itself,  then  noting  that  the  w+  modes  are  much  more  rapidly  evolving  than 
either  the  internal  wave  or  PV  modes  when  c  is  large,  we  may  eliminate  the  acoustic  modes 
by  either  some  form  of  averaging  or  further,  approximate  canonical  transformations  to  the 
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“adiabatic  invariant”  associated  with  the  rapid  motion. 


We  complete  the  cycles  of  canonical  transformation  given  in  this  paper  by  recording  the 
full  expression  of  the  conjugate  coordinates  X(r,  i)  and  7r(r.  t )  in  terms  of  the  normal  modes. 
First  we  recall  that  the  X  and  tt  have  Fourier  expansions 


Xi(r,t)  =  ■  ,  {  Y  <7.(n3,  t)exp(z'Kn3r3) 

sJPqDL  n3;nj.=0 

+  Y  ?,(n,<)exp(f(Ani  ■  r±  + /cn3r3))},  (87) 

a;ni./0 

*.(M)  =  Jjjjzi  Y  P«(n3>0  exp(— t«n3r3) 

V  UL  n3;nj.=0 

+  Y  P.(n,t)exp(-f(An±  •  rj.  +  /cn3r3))}.  (SS) 

n;Uj.*o 


In  terms  of  the  conjugate  variables  >4pv(n),  Bpv( n),  Aj(n),  Bj{ n),  A,i(n),  and  Ba( n)  we 
have  the  total  quadratic  Hamiltonian 


+  +  ul(n)A,(n)A,(~n)  (89) 


—  2j  An±Bpy(n) 


\jul(n)-ujl(n) 


\  ^+(n)  —  wl(n) 


where  we  have  combined  the  two  forms  for  nx  =  0  and  nx  i1  0,  noting  that  at  nx  =  0, 
u>+  =  7^,  and  u ?_  =  /2.  Also  we  have  used  the  notation  that  A,(nx  =  0,n3)  =  a,(n3)  and 
B,(ni  =  0,n3)  =  6i(n3)  to  conform  to  our  earlier  symbols.  The  equations  of  motion  at  linear 
order  are  those  given  in  Eq.(  82). 

Finally  we  record  the  connections  between  the  <7,(n)  and  p,( n)  and  the  normal  mode 
coordinates.  First  for  nx  =  0  we  have 


9i(n3)  =  aPV(n3)- 


bi{n3) 

f 


(90) 
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(91) 


bpyjn 3) 


<?2(n3)  =  a/(n3) - j - , 

PiM  =  |[/a/(n3)  +  W(n3)], 

P2M  =  ^[fa’PV(n3)  +  b]{n3)]. 
For  the  more  complex  case  when  ni  ^  0  we  have  the  following 


?i(“) 

92(n) 

9s(n) 


■(niQjCn)  -  n2Qi(n)), 


=  t — (niQi(n)  +  n3Q2(n)), 

=  <?3(n), 


where 


Qt(  n) 

Q3(n) 

Also  we  have 


^pv(n), 

*"j-  / 

>/w+(n)  -w_(n)  ' 

Anxc(<7!  -  tcr2(n3)) 


-  -v2 

'n3 


^(n) 


A/(n) 


w|(n)-u;2(n)  \Joi2(n)-  723  /)£, -u£(n) 


Pi(n)  = 

tA(«iP2(n)  n2P:(n))  2Anx^1^1^  n2<^2(n)l> 

(100) 

Pa(n)  = 

iA(n1P1(n)  +  n2P2(n))  /  [rijQ^n)  +  n^fn)], 

ZaTIj_ 

(101) 

Pa(n)  = 

^s(n), 

(102) 

and  finally 

^i(n)  =  5pv(n), 

Pj(n)  =  1 — rfh  2 7 '  G/w+<n>  -  ii, 

Xn±Jujl(n)  -  ul(n)  ' 


(103) 
, (104) 
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ft(n)  = 


(wl(n)  -7* ,)(*>« ,  -wi(n)) 


Anxc(<ri  -  icr2(n3))  (n)  -  u>£(n) 


x  (\Ai  - ^(njfl^n)  -  ^'(n)  - ij, «/(»)) 


(105) 


These  connections  summarize  the  set  of  canonical  transformations  we  have  developed  to 
transform  the  Lagrangian  fluid  variables  X(r,  t)  and  ?r(r,  t)  to  normal  modes  of  the  linearized 
problem  of  stratified  flow.  Using  these  expressions  in  the  full  nonlinear  Hamiltonian  Eq.(  83) 
will  allow  the  study  of  the  interaction  among  the  PV,  I,  and  A  modes  to  all  orders  while 
keeping  the  physical  identification  we  have  established.  An  analysis  of  several  of  the  issues 
and  suggestions  made  here  will  be  presented  in  our  subsequent  work. 
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